To Appear in The Astrophysical Journal 

Preprint typeset using L^T^X style cmulatcapj v. 6/22/04 



in 
o 
o 

(N 

C 
3 



(N 

> 



in 

O 

in 
o 

Oh! 

6 



13 



THE SPECTRAL EVOLUTION OF TRANSIENT ANOMALOUS X-RAY PULSAR XTE J1810-197 

E. V. GOTTHELF AND J. P. HALPERN 
Columbia Astrophysics Laboratory, Columbia University, 550 West 120 th Street, 
New York, NY 10027-6601; eric@astro.columbia.edu 
To Appear in The Astrophysical Journal 

ABSTRACT 

We present a multi-epoch spectral study of the Transient Anomalous X-ray Pulsar XTE J1810— 197 
obtained with the Newton X-Ray Multi-Mirror Mission (XMM-Newton) . Four observations taken over 
the course of a year reveal strong spectral evolution as the source fades from outburst. The origin of 
this is traced to the individual decay rates of the pulsar's spectral components. A two-temperature 
fit at each epoch requires that the temperatures remains nearly constant at kT\ = 0.25 keV and 
kT-2 = 0.67 keV while the luminosities of these components decrease exponentially with t\ = 900 days 
and T2 = 300 days, respectively. The integrated outburst energy is E\ = 1.3 x 10 42 d\ 5 kpc ergs and 

E 2 = 3.9 x 10 42 d\ 5 kpc ergs for the two spectral components, respectively. One possible interpretation 
of the XMM-Newton observations is that the slowly decaying cooler component is the radiation from a 
deep heating event that affected a large fraction of the crust, while the hotter component is powered by 
external surface heating at the foot-points of twisted magnetic field lines, by magnetospheric currents 
that are decaying more rapidly. The energy-dependent pulse profile of XTE J1810— 197 is well modeled 
at all epochs by the sum of a broad pulse that dominates in the soft X-rays and a narrower one at 
higher energies. These profiles peak at the same phase, suggesting a concentric emission geometry 
on the neutron star surface. The spectral and pulse evolution together argue against the presence of 
a significant "power-law" contribution to the X-ray spectrum below 8 keV. The extrapolated flux is 
projected to return to the historic quiescent level, characterized by an even cooler blackbody spectrum, 
by the year 2007. 

Subject headings: pulsars: general — stars: individual (XTE J1810-197) — stars: neutron — X-rays: 
stars 



1. INTRODUCTION 

The bright 5.54 s X-ray pulsar XTE J1810-197 is the 
second example of a Transient Anomalous X-ray Pul- 
sar (TAXP) and the first one confirmed by measuring 
a rapid spin-down rate. All of its observed and derived 
physical parameters are consistent with classification as 
an Anomalous X-ray Pulsar (AXP), one that had an im- 
pulsive outburst sometime between 2002 November and 
2003 January, when it was discovered serendipitously by 
llbrahim et alJ l)2004|) using the Rossi X-ray Timing Ex- 
plorer (RXTE). Its flux was observed to be declining with 
an exponential time constant of 269±25 days from a max- 
imum of F x (2 - 10 keV) w 6xl0~ n ergs cm" 2 s" 1 . The 
source was then localized precisely using two Target of 
Opportunity (T oO) observations with the Chandra X-ray 
Obs ervatory by Gotthelf et alJ (|2004t hereafter Paper I) 
and llsrael et~al. l|2004f l. In comparison, archival detec 
tions by several X-ray satellites indicate a long-lived qui 
escent baseline flux of F x (0.5 - 10 keV) « 7 x 10 
ergs cm~ 2 s _1 lasting at least 13 years and possibly 
for 23 years prior to the outburst (Paper I). Fading 
of an IR source within the Chandra error circle, simi- 
lar to ones associated with other AXPs, confirmed it s 
identification with XTE J1810-197 l|Rea et all2004albl h 
The X-ray spectra and pulse profiles from three observa- 
tions obtained with XMM-N ewton during the dec l ine of 
the outburst were studied bv lHalpern fc Gotthelfl (2005, 
hereafter Paper II). The short duty cycle of activity of 
XTE J1810— 197 suggests the existence of a significant 
population of as-yet unrecognized, although not neces- 
sarily undetected, young neutron stars (NSs). 
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In this paper we present a new XMM-Newton observa- 
tion of XTE J1810-197. This data set, combined with 
previous XMM-Newton observations acquired over the 
last year, allows us to characterize the spectral evolution 
of a Transient AXP. We show that the decay rates of the 
two fitted X-ray spectral components are distinct, and 
can be extrapolated in time to approach the previous 
quiescent spectrum. For consistency with Papers I and 
II, we express results in terms of a maximum distance 
of d = 5 kpc to the pulsar. However, there is reason to 
believe that XTE J1810— 197 is significantly closer than 
this, and we take into account a more realistic estimate 
of 2.5 kpc in our discussion of proposed physical models 
for the X-ray emission and outburst mechanism. 

2. OBSERVATIONS 

A fourth XMM-Newton observation of 
XTE J1810-197 was obtained on 2004 September 
18. The previous three observations are described in 
Paper II. We use the data collecte d with the Eu r opean 
Photon Imaging Camera (EPIC, iTurner et al.l 12003(1 
which consist of three CCD imagers, the EPIC pn 
and the two EPIC MOSs, each sensitive to X-rays in 
the 0.1 — 12 keV energy range. In the following we 
concentrate on data taken with the EPIC pn detector, 
which provided a timing resolution of 48 ms in "large 
window" mode, for ease of comparison with the earlier 
data sets. The fast readout of this instrument ensures 
that its spectrum is not affected by photon pileup. Data 
collected with the two EPIC MOS sensors used the 
"small window" mode and "timing" mode, providing a 
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Fig. 1. — Period evolution of XTE J1810— 197 using timing mea- 
surements obtained with XMM-Newton as listed in Table HI The 
dashed line is a chi-square fit yielding the indicated mean value of 
P. Error bars are 95% confidence and indicate probable deviations 
from a constant spin-down rate. 



time resolution of 0.3 s and 1.5 ms, respectively. 

For the new data set, we followed the reduction and 
analysis procedures used for the previous XMM-Newton 
observations of XTE J1810— 197, as outlined in Paper 
II. The new observation is mostly uncontaminated by 
flare events and the filtered data set resulted in a total 
of 26.5 ks of good EPIC pn exposure time (24.4 ks live 
time). We checked for timing anomalies that were evi- 
dent in some previous EPIC pn data sets and found none. 
Photon arrival times were converted to the solar system 
barycenter using the Chandra derived source coordinates 
R.A. 18 h 09 m 51. s 08, deck -19°43'51."7 (J2000.0) given in 
Paper I. 

2.1. Spin-down Evolution of XTE J1810-197 

The barycentric pulse period of XTE J1810— 197 mea- 
sured at each XMM-Newton observing epoch is given in 
Tabled These are derived from photons obtained with 
both the EPIC pn and MOS cameras, with the exception 
of the short 2003 October 12 observation for which only 
EPIC pn data of sufficient time resolution is available. 
The errors are the 95% confidence level determined from 
the Z\ test. As shown in Figure^ the four period mea- 
surements can be fitted to yield a mean spin-down rate of 
P = (8.1 ±0.7) x 10~ 12 s s^ 1 over the year-long interval. 
This implies a characteristic age r c « 10,800 yr, surface 
magnetic field B s w 2.1 x 10 14 G, and spin-down power 
E » 1.9 x 10 33 erg s -1 , comparable to the earlier val- 
ues l)Ibrahim et alJl2004 and Paper II). Deviations from 
a cons tant P are evident, however, since llbrahim et akl 
(|2004[) fitted values in the range (1.1 - 2.2) x 10" 11 s s _1 
in the first 9 months of the outburst. 

As discussed in paper II (§2.2), timing glitches 
in AXPs can result in large increases in their pe- 
riod derivatives \AP/P ~ 1; e.g. , I E 2259+586 
llwasawa. Kovama. fc Halnernl l)1992tk iKasoi et al. 
(120031) and 1RXS J170849.0-400910 lDaH'Osso et al 
( 2003)]. Accordingly, it is possible that XTE J1810-197 
experienced a glitch and an increase in P at the time 
of its outburst, and that P is now relaxing to its long- 



term value. But since there is no prior ephemeris for 
XTE J1810— 197 in its quiescent state, we do not know if 
its outburst was triggered by a glitch. Alternatively, the 
spin-down torque could have been enhanced in the early 
stages of the outburst by an increase over the dipolc 
value of t he magnetic field strength at the speed- of-light 
cylinder (|Thompson. Lvutikov. fcKu lkarni 2 002[). or bv 
a particle wind and Alfven waves (iThonrosmTfe Blaesl 
119981 lHarding. Contonoulos. & Kazana^Tl999y . effects 
that are expected to decline after the first few months. 

2.2. Spectral Analysis and Results 

XMM-Newton observations of XTE J1810-197 have 
shown that its spectrum is equally well fitted by a power- 
law plus blackbody model, as commonly quoted for 
AXPs, or a two-temperature blackbody model. In Pa- 
per II we argued that the two-temperature model is more 
physically motivated, while the power-law plus black- 
body model suffers from physical inconsistencies. As 
we shall show, the new data bolster these arguments, so 
we concentrate mainly on the double blackbody model 
in this work. Table ^ presents a summary of spec- 
tral results from all four XMM-Newton observations of 
XTE J1810-197. 

As with the earlier data, the 2004 September source 
spectrum was accumulated in a 45" radius aperture 
which encloses > 95% of the energy. Background 
was taken from a circle of the same size displaced 2.'3 
along the readout direction. The spectra were grouped 
into bins containing a minimum of 400 counts (includ- 
ing background) and fitted using the XSPEC package. 
For this analysis the column density is held fixed at 
Nh = 6.5 x 10 21 cm -2 , the value determined from pre- 
vious fits, which are all consistent. The best fit to the 
two-temperature blackbody model yields temperatures 
of kTx = 0.25 ± 0.01 keV and kT 2 = 0.67 ± 0.01 keV 
with a fit statistic of xt — 1-2 for 188 degrees of freedom 
(Fig.EJ). Each of these two temperatures, which we refer 
to as "warm" and "hot", respectively, remained essen- 
tially the same, nominally to within 4% of that reported 
for the first XMM-Newton observation (see Table EJ. 
However, the warm and hot blackbody luminosities de- 
clined by 33% and 70%, respectively, in one year. With 
the temperature of the hot component remaining essen- 
tially constant, its decline in luminosity is attributable 
to a decrease in its emitting area, A 2 . The most recent 
value, A 2 = 2.1 x 10 11 a% cm 2 , is « 1% of the NS surface 
area. In the case of the warm component, the decay in 
flux is modest, so that it is not yet possible to decide 
within the errors whether the temperature or the area 
is the primary variable. The most recent area measure- 
ment of the warm component, A\ = 9.1 x 10 12 d\ cm 2 , is 
w 50% of the NS surface area. 

Phase-resolved spectroscopy using the new XMM- 
Newton data set is compared here to that reported in 
Paper II. As before, we fit for the intensity normalization 
in each of 10 phase bins with the two temperatures and 
column density fixed at the phase-averaged value listed 
in Table ^ This is equivalent to determining the rela- 
tive projected area of emission as a function of rotation 
phase. The fits are again found to be each statistically 
acceptable as a set, so no significant test can be made 
for variations in additional parameters such as the tem- 
peratures or column density. The results of the phase- 
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ENERGY (KEY) ENERGY (KEV) 

Fig. 2. — Comparison of XMM-Newton EPIC pn spectra of XTE J1810-197 obtained at (a) the earliest epoch (2003 September) and 
(b) the latest epoch (2004 September), fitted with a two- temperature blackbody model as described in the text and specified in Table HI 
Although the temperatures of the blackbody components have not changed between the two epochs, the flux of the hot component (BB2) 
has decayed faster than the warm one (BB1). Also shown are the residuals from the best-fit models. 



TABLE 1 

XMM-Newton Spectral and Timing Results for XTE J1810-197 



Parameter 




2003 Sep 8 


2003 Oct 12 


2004 Mar 11 


2004 Sep 18 


N H (10 22 cm- 2 ) a 




0.65 ±0.04 


0.65 ±0.04 


0.65 ±0.04 


0.65 (fixed) 


fcTi (keV) 




0.26 ± 0.02 


0.29 ±0.04 


0.27 ±0.02 


0.25 ±0.01 


kT 2 (keV) 




0.68 ±0.02 


0.71 ±0.03 


0.70 ±0.01 


0.67 ±0.01 


Area Ai (cm 2 ) 




1.1 x 10 13 


6.6 x 10 12 


6.8 x 10 12 


9.1 x 10 12 


Area A2 (cm 2 ) 




6.4 x 10 11 


5.1 x 10 11 


2.9 x 10 11 


2.1 x 10 11 


BB1 Flux b 




4.2 x 10" 12 


5.4 x 10~ 12 


3.5 x 10~ 12 


2.6 x 10" 12 


BB2 Flux b 




3.5 x 10" 11 


3.0 x 10- 11 


1.8 x 10- 11 


1.0 x 10- 11 


Total Flux b 




3.93 x 10- 11 


3.84 x 10" 11 


2.13 x lO" 11 


1.29 x 10" 11 


^BBi(bol) (ergs s" 


1)C 


5.2 x 10 34 


5.1 x 10 34 


3.9 x 10 34 


3.5 x 10 34 


^BB2(bol) (ergs s" 


1)C 


1.4 x 10 35 


1.3 x 10 35 


7.2 x 10 34 


4.2 x 10 34 


xS(dof) 




1.1(187) 


1.1(84) 


1.1(194) 


1.2(188) 



EPIC pn exposure (ks) 11.5 6.9 17.0 26.5 

EPIC pn live time (ks) 8.1 6.2 15.8 24.4 

Off-axis angle (arcmin) 0.0 8.8 0.0 0.0 

Count rate (s" 1 ) 11 10.6 4.8 5.8 3.4 

Epoch (MJD/TDB) C 52890.5642044 52924.0000320 53075.4999960 53266.4999776 

Period (s) f 5.539344(19) 5.53943(8) 5.539425(16) 5.539597(12) 



Note. — Uncertainties in spectral parameters are 90% confidence for two interesting 
parameters. 

a Parameter derived from a linked fit to all epochs. 
b Absorbed 0.5-10 keV flux in units of ergs cm -2 s — x . 
c Bolometric luminosity assuming d = 5 kpc. 

d Background subtracted EPIC pn count rate corrected for detector dead-time. 
c Epoch of phase zero in Figure 151 

includes EPIC MOS data where available. 95% confidence uncertainty in parentheses. 



resolved spectral fits are shown in Figure Compari- 
son between the data sets taken a year apart shows that 
the modulation with phase of each blackbody component 
has remained steady to within the measurement errors. 
The phase alignment of the two temperature components 
has remained the same, and the pulses peak at the same 
phase. This is consistent with the picture of a small, hot 
region surrounded by a warm, concentric annulus that 
occupies ~ 1/2 the surface area of the star. Neither com- 
ponent disappears at any rotation phase. In particular, 
if the hot component were completely eclipsed, the spec- 
tral decomposition in Figure [21 indicates that the light 
curves at E > 3 keV should dip to zero; clearly they do 



not. 

2.3. Long-term Flux Decay 

With the set of XMM-Newton measurements spanning 
a year, a more accurate characterization of the X-ray flux 
decay of XTE J1810-T97 is possible than that rep orted 
in the initial RXTE study of llbrahim et all {2004). The 
bolometric flux over time, shown in Figure 0] reveals 
a new and complex behavior. While the RXTE data 
were consistent with an exponential decay of time con- 
stant r = 269 ± 25 days, the subsequent XMM-Newton 
measurements show that the inferred bolometric lumi- 
nosity of the two spectral components are not declin- 
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Fig. 3. — Contribution of the two blackbody components as a 
function of rotation phase at two epochs. ( Top) 2003 September. 
(Bottom) 2004 September. The temperatures kT\ and kT2 are 
held fixed at the values listed in Table while the normalization 
constants are fitted. Inset: For each epoch, the ratio of blackbody 
normalization constants as a function of rotation phase. 



ing at the same rate. From the data presented in Ta- 
ble ^ we find that the luminosity of the hotter tempera- 
ture component falls exponentially with r 2 = 300 days, 
while the warm component decreases with t\ = 900 days, 
which, although less reliably characterized, is clearly 
longer than t 2 . The shorter time constant is very close 
to the one describing the RXTE data. Furthermore, the 
RXTE spectrum itself is fitted bes t by a single tempera - 
ture blackbody of kT = 0.7 keV (jRoberts et al.ll200^) . 
without a power-law contribution, which is consistent 
with our interpretation of the XMM-Newton spectrum. 
Since RXTE is sensitive only at energies > 2 keV where 
the hotter blackbody component dominates, the over- 
all X-ray decay from the beginning of the outburst to 
the present appears consistent with separate exponen- 
tial time constants corresponding to the two distinct 
thermal-spectrum components. We also note that it is 
not possible to fit an alternative power-law temporal de- 
cay to the hot blackbody flux; such a model would require 
a decay index that steepens with time. 

As shown in Figure we expect that the X-ray flux of 
XTE J1810— 197 will return to its historic quiescent level 
by the year 2007. If the quiescent spectrum is to match 



that observed before the current outburst, then it would 
likely resemble the ROSAT observation of 1992 March 
7. Although of relatively poor quality, the ROSAT spec- 
trum can be reasonably well fitted with a single black- 
body of kT = 0.18 ± 0.02 keV covering 1.2 x 10 13 d\ cm 2 
and L BB (bol) = 1.3 x 10 34 rf| ergs s^ 1 (Paper I).' This 
blackbody is significantly cooler and larger than the out- 
burst warm component of kT\ = 0.25 keV and A\ = 
9.1 x 10 12 <ig cm 2 , so it should be present even now, 
although masked by the fading outburst emission. We 
estimate that the kT = 0.18 keV blackbody spectrum 
measured with ROSAT would contribute 30% of the lat- 
est (2004 September) XMM-Newton measured flux at 
0.5 keV. Therefore, we expect that such a cooler com- 
ponent will begin to dominate the soft X-ray spectrum 
in late 2006, consistent with the projection of the flux in 
Figure 0] 



2.4. Pulse Shape Evolution 

The energy-resolved folded light curves from the 2003 
September 8 XMM-Newton observation (Fig. [5^) show 
that the pulse peak, in general, is somewhat narrower 
than a sinusoid, an effect that is more pronounced at 
higher energy. Examination of the light curves measured 
one year later (2004 September 18; Fig.[5ji) shows a clear 
energy-dependent change. This is easily seen by forming 
the ratio of the curves taken at the two epochs, normal- 
ized by their count-weighted mean. These ratio curves, 
presented in Fig. [fJl show changes in amplitude of up to 
20%. Specifically, the pulse shape of the lower-energy 
light curves has evolved to a smaller amplitude, more si- 
nusoidal profile. This change is highly significant as de- 
termined by the reduced xt (24 degrees-of-freedom) for 
the null hypothesis of a flat line (see labels on Fig. [fjj). 
In particular, the most pronounced variation is found in 
the 1.0-1.5 keV band with xt = 6.3 corresponding to 
an infinitesimal probability of the pulse shape being con- 
stant. In contrast, the pulse shape of the higher energy 
bins remained essentially unchanged, as evidence by their 
xt values near unity. This trend is reflected in the two 
intervening observations (Fig. [3J; and[S|i) as well. 

The pulse profile must be made up of some combina- 
tion of pulsed emission intrinsic to each of the two mea- 
sured spectral components. In any given energy band, 
the changing relative contribution of the soft and hard 
spectral components forces the combined pulse shapes 
to evolve in time. One way of quantifying this change 
is to define the pulsed fraction f p , the fraction of flux 
above the minimum in the light curve. This is found 
to increase smoothly with energy from 34% at less than 
1 keV to « 51% above 5 keV at the earlier epoch (see 
Tableland Fig.[SJ). At the latest epoch reported herein, 
the pulsed fraction in the lower energy band decreased 
to 26% while it remained essentially unchanged at the 
higher energy. 

The simplest hypothesis for modeling the light curve 
evolution is to assume that the intrinsic light curve of 
each spectral component is steady in time. If so, the 
pulsed fraction of the light curve at a given energy and 
epoch f p {E,t) can be expressed as a linear combination 
of two light curves of fixed pulsed fractions fsBi and 
/b B2 weighted by the ratio of counts from the respective 
spectral component NsBi{E,t) and NsB2(E,t), 
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Fig. 4. — The bolometric luminosity and pulsed fraction of 
XTE J1810-197 as a function of time. (Top) The bolometric 
luminosity of the individual components of the two-temperature 
blackbody model (crosses) derived from spectral fits to XMM- 
Newton data with the temperatures fixed (see text). The fitted 
c-folding times are n = 900 d and T2 = 300 d for the warm and 
hot blackbody components, respectively. The fitted quantities have 
been extrapolated to the (lc) quiescent range measured in Paper 1 
(cross-hatched area). (Bottom) Pulsed fraction measurements in 
the 1.0 — 1.5 keV energy band. The solid line represents the pulsed 
fraction modeled in §2.4 from the fitted decay of the blackbody 
components shown in the top panel. 



the double blackbody spectral model tabulated in Ta- 
bleau Treating the ratio R{E, t) as the independent vari- 
able, the best fit yields Jbb2 = 52 ± 3% for the hotter 
component and Jbbi = 26 ± 2% for the cooler one. The 
fit statistic is x2 = 1-17 for 22 degrees-of- freedom, corre- 
sponding to a probability p(> 1.17) = 0.27. The lower 
panel of Fig. 21 shows the best fit model for the decay of 
the the 1.0-1.5 keV pulsed fraction over time and its ex- 
trapolation to an epoch dominated by the warm spectral 
component. On the other hand, if we apply the black- 
body plus power-law model, the corresponding values of 
R(E,t), also listed in Table are completely different 
and the result is xt — 1-70 for 22 degrees-of-freedom, 
corresponding to a probability p(> 1.70) = 0.021, 13 
times less likely relative to the double blackbody model, 
for the assumption of fixed intrinsic pulsed fractions. 

Evidently the pulsed fraction of the hotter component 
is much higher than for the cooler component. This 
difference would account for the gradual shift from a 
sharper, more triangular pulse shape found at higher en- 
ergies, to a rounder and more symmetric sinusoidal light 
curve seen at lower energies. As was hypothesized above, 
the decrease in pulsed fraction with time at low energies 
follows from the fact that the hot spectral component 
is decaying more rapidly than the cooler one, and so its 
contribution to the pulsed fraction in the low-energy bin 
is declining. To the extent that we are not able to dis- 
cern differences greater than « 4%, which represents the 
measurement uncertainty, the pulsed fractions Jbbi and 
Jbbi are individually independent of energy. However, 
we also cannot rule out a model in which Jbbi itself 
increases with increasing energy. 

It is useful to model the shape of the light curves of 
XTE J1810-197 in detail to further quantify their evo- 
lution. Overall, the pulse profiles suggest some linear 
combination of sinusoidal and triangular functions, for 
the soft and hard X-rays, respectively. In the following, 
we model the photon counts as a function of phase, N(<fi), 
at a given energy E and epoch t by the two-component 
model, 



f P (E,t) 



_ Sb 



N B Bi(E,t)+f BB2 N BB2 (E,t) 
N BB1 (E,t)+N BB2 (E.t) 



= fBBl[l- R(E,t)}+ f B B2R(E,t), (1) 

where R(E,t) = N BB2 {E,t)/[N BB1 {E,t) + N BB2 {E,t)] 
is the normalized flux ratio. This expression is only valid 
in the case that the minima of the two light curves coin- 
cide, as is evident for the XTE J1810-197 profiles. For 
an exponential decay of the component luminosities, the 
time-dependent model for the pulsed fraction is simply, 



where 



N((t>- E, t) = Nsifc E, t) + N T (4>- E, t), 



N s (<t>; E, t) = a(E, t) [ 1 + cos(<£ -&?)] + 7s(E, t) 



and 



N T {<t>;E,t) 



p(E,t)[\-2\(h- 
1t(E, t), 



cj> T \/5(E,t)]+ 7T (E,t), 
if\<f>-<fr\<8/2 
if\<t>-<fa\>5/2 
(3) 



f P (E,t) 



fBBi + .fBB2r{E,t )e- ( - t - t ^/ T " 
l + r(£;,t )e-(*- t °)/^ 



(2) 



where r(E,t a ) = N B m(E,t )/N BB1 (E,t ) is the flux 
ratio at time t and Td = (1/t2 — l/ri) -1 = 450 d is the 
differential time constant. 

We can test this hypothesis directly by fitting for Jbbi 
and Jbbi using a joint least-squares fit to the set of 24 
(six energy bands at four epochs) measured pulsed frac- 
tions f p (E,t) given in Table Hand the ratio R(E,t) for 



Here, a is the amplitude of the pulsed signal and 7s 
represent the "unpulsed," or minimum level, for the si- 
nusoidal component, (3 and 7t are the analogous param- 
eters for the triangular component, and S is the duty 
cycle (full width) of the triangular pulse. 

Initial fits to the light curves indicate that the model 
can be constrained by fixing the triangular width S to a 
single value for all observations. The same proves true 
for the relative phase 4>s — 4>t between the peak fluxes of 
the two components. Accordingly, we used a bootstrap 
approach to determine these shape parameters. After 
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Fig. 5. — Energy-dependent pulse profiles of XTE J1810— 197 obtained with the XMM-Newton EPIC pn detector. This page shows 
profile from epochs (a) 2003 September 8 and (6) 2003 October 12. Also shown is the best fit to the two-component model for the pulse 
profile (solid line) described in the text (see §2.4). The individual contibutions are plotted for the triangular component (Phase — 1) and 
the sinusoidal component (Phase 1 — 2). 
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Fig. 5. Continued - Energy-dependent pulse profiles of XTE J1810— 197 obtained with the XMM-Newton EPIC pn detector for (c) 2004 
March 11 and (d) 2004 September 18. The epochs of phase zero are given in Table^ Background has been subtracted. Pulsed fraction at 
low X-ray energies has decreased in time, while it has remained essentially unchanged at high energy. 
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TABLE 2 

XMM-Newton Pulse Profile Fit Results for. XTE J1810-197 



Energy 


2003 




2004 




Range 


Sep 8 


Oct 12 


Mar 11 


Sep 18 


(keV) 


(%) 


(%) 


(%) 


(%) 



Pulsed Fraction f p and (Fit Statistic x?) a 



0.5 - 1.0 34 ±4 (0.9) 33 ±6 (1.2) 31 ± 4 (0.5) 26 ± 3 (0.9) 

1.0 - 1.5 44 ±2 (0.6) 43 ±3 (1.5) 40 ± 2 (0.8) 36 ± 2 (1.2) 

1.5 - 2.0 48 ±2 (1.5) 48 ±3 (1.3) 46 ± 2 (1.5) 42 ± 2 (1.6) 

2.0-3.0 51 ±2 (1.4) 51 ±3 (1.2) 50 ± 2 (0.9) 45 ± 2 (1.9) 

3.0- 5.0 54 ±2 (1.2) 55 ± 3 (1.0) 50 ± 2 (1.0) 50 ± 2 (1.5) 

5.0-8.0 51 ±6 (2.0) 62 ±5 (0.9) 48 ± 6 (1.5) 47 ±6 (1.6) 



Normalized Sinusoidal Component Amplitude a 



0.5 - 


1.0 


28 ±4 


27 it 7 


30 ±4 


22 ±3 


1.0 - 


1.5 


40 ±2 


37 ±3 


37 ±2 


31 ±2 


1.5 - 


2.0 


39 ±2 


44 ±3 


41 ±2 


39 ±2 


2.0 - 


3.0 


40 ±2 


39 ± 3 


42 ±2 


33 ±2 


3.0 - 


5.0 


39 ±3 


37 ±4 


34 ±3 


33 ±3 


5.0- 


8.0 


35 ±8 


34 ±8 


15 ±8 


06 ±8 



Normalized Triangular Component Amplitude /35/2 



0.5 - 


1.0 


6±3 


6±4 


1±2 


4±7 


1.0 - 


1.5 


5±2 


6±2 


4± 1 


4±6 


1.5 - 


2.0 


10 ±2 


4±2 


5±1 


4±6 


2.0 - 


3.0 


11 ±2 


11 ±3 


8±1 


11 ±6 


3.0 - 


5.0 


14 ±2 


18 ±3 


16 ±2 


17 ±5 


5.0 - 


8.0 


14 ±5 


27 ±6 


32 ± 5 


40 ±5 



Normalized Unpulsed Amplitude 7 



0.5 - 


1.0 


67 ±4 


66 ±3 


69 ±3 


74 ± 2 


1.0 - 


1.5 


57 ±2 


56 ±2 


59 ±2 


64 ±1 


1.5 - 


2.0 


52 ±2 


52 ± 2 


54 ±2 


58 ± 2 


2.0 - 


3.0 


49 ±2 


49 ±2 


50 ±2 


55 ± 2 


3.0 - 


5.0 


45 ±2 


46 ±2 


49 ±2 


50 ± 2 


5.0- 


8.0 


48 ±5 


38 ±4 


52 ± 5 


52 ±4 



Note. — Results of fits to the two-component sinusoidal plus 
triangular model for the pulse profile described in §2.4. The fixed 
model parameters are (f>s — 4>T = 0.0, 5/27T = 0.53 cycle, and at each 
epoch, the absolute phase of the sinusoidal function. 

a Uncertainties in the pulsed fractions are at the 1-cr (68% confi- 
dence level) for three interesting parameters (a, (3, 7). 



finding a global fit by iteration, the triangular width is 
set at S/2tt = 0.53 cycles. We then determined (j>s — 4>t, 
which is found to be consistent with zero. Finally, the 
absolute phase 0s at each epoch was fixed. We therefore 
conclude that there is no need to allow energy or time 
dependence for the individual component pulse shapes. 

With the shape of the pulse components fixed, we fit 
the restricted model with just a, /?, and 7 (= 75 + 7t, 
since the 7's are indistinguishable in this fit) as free pa- 
rameters. This model fitted to all 24 individual pulse 
profiles yields a reduced xt ranging from 0.6 — 2.0 with a 
typical value of xt — 1-2, corresponding to a probability 
°f Xv) — 0-23 for 22 degrees-of-freedom. The best- 
fit model parameters, derived pulsed fractions, and fit 
statistic for each light curve are presented in Table |2 A 
clear trend is seen as the pulsed fraction decreases with 
time. The best-fit models are shown in Figure overlaid 
on plots of the energy-resolved pulsed profiles for the 
four epochs. We checked that a single component model, 
either a pure sinusoidal or triangular pulse, does not ad- 
equately characterize the data based on an unacceptably 
large xt f° r many of the fits. 



TABLE 3 

Component Flux Fraction for XTE J1810— 197 



Energy 2003 2004 



Rang 


;e (keV) 


Sep 8 


Oct 12 


Mar 11 


Sep 18 




R(E,t) 


= N BB2 


/(Nbbi ' 


fW B B2) a 




0.5 - 


1.0 


0.389 


0.332 


0.285 


0.218 


1.0 - 


1.5 


0.591 


0.499 


0.466 


0.399 


1.5 - 


2.0 


0.800 


0.702 


0.696 


0.661 


2.0 - 


3.0 


0.941 


0.886 


0.896 


0.893 


3.0 - 


5.0 


0.996 


0.988 


0.991 


0.992 


5.0 - 


8.0 


1.000 


1.000 


1.000 


1.000 



R(E,t) = N PL /(N PL + N BB ) b 



0.5 - 


1.0 


0.928 


0.911 


0.925 


0.952 


1.0 - 


1.5 


0.755 


0.718 


0.753 


0.807 


1.5 - 


2.0 


0.533 


0.480 


0.524 


0.567 


2.0 - 


3.0 


0.347 


0.291 


0.331 


0.340 


3.0 - 


5.0 


0.247 


0.191 


0.223 


0.200 


5.0 - 


8.0 


0.376 


0.278 


0.325 


0.261 



*R(E, t) = N BB2 1 (N BB1 + N BB2 ) is the fraction of 
counts in a given energy band from the hotter com- 
ponent of the two-temperature blackbody spectrum 
based on the spectral model presented in Table H\ 

b R(E,t) = N PL /(N PL + N BB ) is the fraction of 
counts in a given energy band from the power-law 
component of the best fit power-law plus blackbody 
model blackbody spectrum (not shown). 



It is possible to further reduce the number of free pa- 
rameters in the fit by constraining the coefficients a, (3, 
and 7 so that the phase-averaged photon count ratios 
Ns(E,t)/NT(E,t) are identical with the ratios derived 
from the two-temperature blackbody spectral fits in each 
energy interval at each epoch. That is, we test the 
hypothesis that the warm blackbody is responsible for 
only the sinusoidal pulse and the hot blackbody for only 
the triangular pulse by requiring Ns(E,t)/NT(E,t) = 
N Bm {E,t)/N B B2{E,t) at all epochs. We find that, 
while reasonable fits to the light curves are possible 
for the lower and highest energy bands, fits in the 
1.5 — 5.0 keV range proved unacceptable. This argues 
against a simple one-to-one correspondence between the 
spectral and the chosen light curve component decompo- 
sition. It is possible that we have not started with the 
correct spectral or light curve basis pairs or that the basic 
pulse profiles are each energy dependent. It is also possi- 
ble that there may be a third, unfitted softer component 
associated with the quiescent flux, as seen by ROSAT 
(see §2.3). Continuing observation of the spectrum as 
the source fades should clarify the relationship between 
spectral and light curve components. 

3. INTERPRETATION 

3.1. Spectral Evolution: Constraining Models 

While the 0.5 — 8 keV X-ray spectra of AXPs (including 
XTE J1810— 197) clearly require a two-component fit, it 
is not possible to prove from spectra alone that either 
of two very different models, namely a power-law plus 
blackbody, or a two-temperature blackbody, is correct. 
For XTE J1810— 197, either model yields an acceptable 
chi-square when fitted to the X-ray spectrum. However, 
in Paper II we presented physical arguments against the 
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properties of the particular power law that results from 
the power-law plus blackbody fit. Unlike some rotation- 
powered pulsars whose spectra are fitted by hard power 
laws plus soft blackbodies, in AXPs the roles of these two 
components are reversed. For XTE J1810— 197, most of 
the X-rays belonging to the T ~ 3.7 power-law compo- 
nent are lower in energy than those fitted by the black- 
body (see Figure 5a of Paper II). This steep power law 
must turn down sharply just below the X-ray band in 
order not to exceed the faint, unrelated IR fluxes. As 
shown in Paper II, this would be difficult to achieve in 
a synchrotron model. And in a Comptonization model, 
where are the seed photons? Even if such a sharp cut- 
off existed at the low end of the EPIC energy band, 
a model in which the power-law component is due to 
Compton upscattering of thermal X-rays from the sur- 
face does not explain why most of the power-law photons 
have lower energy than the observed blackbody photons 
unless there is a larger source of seed photons at ener- 
gies below 0.5 keV. However, such a seed source would 
have to be thermal emission from the neutron star. Since 
the two-temperature model uses a large fraction of the 
surface area for the 0.5 — 2 keV photons, we prefer to 
regard the X-rays in this band as plausible seeds for, 
rather than as the product of, inverse Compton scatter- 
ing. Supporting this assumption in the specific case of 
XTE J1810-197 are the shapes of its pulse profiles. The 
smooth increase in pulsed fraction with energy, and the 
strict phase alignment at all energies, are unnatural un- 
der the hypothesis of two different emission mechanisms 
and locations. 

The evolving spectral shape and pulse profiles of 
XTE J1810— 197 during its decline from outburst offer 
new and complementary evidence concerning the appro- 
priate spectral decomposition. We find that the the 
pulsed fraction is declining with time at low energies, 
while remaining essentially constant at high energies. 
This fact is consistent with the changing contributions of 
fluxes in the two-temperature fit as the two components 
decline at different rates. The warm and hot compo- 
nents contribute significantly to the soft X-rays, whereas 
the hard X-rays are supplied entirely by the hot compo- 
nent. This is quantified in Table where it is seen that 
the hot blackbody component accounts for 30-60% of the 
flux below 1.5 keV, and > 99% of the flux above 3 keV at 
all times. Under this spectral decomposition, the pulsed 
fraction at energies above 3 keV is large and unchanging 
because one and only one spectral component is present 
at those energies, and it has the larger pulsed fraction 
of the two spectral components. In contrast, the pulsed 
fraction at low energies (< 1.5 keV) is understood to de- 
crease in time because the hot component makes only 
a fractional contribution to its light curve, and its frac- 
tion decreases faster than the warm spectral component, 
which has the intrinsically smaller pulsed fraction. This 
effect is also illustrated in Figure 0] where the modeled 
pulse fraction is compared to the data. 

In contrast, consider how the spectral components 
would contribute to the various energy bands in the 
power-law plus blackbody spectral model. Table 3 shows 
that, unlike the two blackbody model, a fitted power-law 
component accounts for 79% of the soft (< 1.5 keV) X- 
rays at all epochs, while decreasing its contribution to 
the hard (> 3 keV) X-rays over time. In this model, the 
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Fig. 6. — Comparison of the XMM-Newton EPIC pn pulse 
profiles of XTE J1810— 197 obtained one year apart in six en- 
ergy bands. Displayed is the ratio in each band between the 2003 
September 8 and 2004 September 18 profiles shown in Fig. 5, nor- 
malized by the count-weighted mean. The profiles show more sys- 
tematic change at lower energies as indicated by the ratio of the 
amplitudes (vertical scale). 
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pulsed fraction should remain constant at low energies 
because only the power law is present there, while at high 
energies the pulsed fraction should vary because both the 
blackbody and the power-law components are present, in 
varying proportions. This spectral decomposition would 
drive an evolution of the pulse shapes that is opposite of 
what is observed. Thus, we find that the detailed evo- 
lution of the X-ray emission from XTE J1810- 197 fur- 
ther supports the assumption of a purely thermal surface 
emission model, and leaves no evidence of a steep power 
law in the 0.5 — 8 keV band, as is commonly fitted to 
individual observations of AXPs. 

Our decomposition of the pulsed light curves is also 
consistent with the absence of modulation in the ROSAT 
observation of 1992 March 7, to an observed upper limit 
of 24% in the 0.2-2.0 keV PSPC energy band. Since the 
ROSAT spectrum is fitted by an even cooler blackbody 
of kT — 0.18 keV, it can be expected that the pulsed 
fraction in quiescence will fall below the 25% value fitted 
to the warm blackbody in XMM-Newton data once that 
component fades away. 

3.2. Re-evaluation of the Distance to XTE 31810—191 

For consistency with Papers I and II, we have parame- 
terized results here in terms of a distance of 5 kpc. How- 
ever, this was considered an upper limit based on an 
H I absorption kinematic dist ance to the neighbo ring su- 
pernova remnant Gil. 2— 0.3 l)Oreen et al.lll98S^ and its 
X-ray measured colu mn density TVh ~ 1.4 x 10 22 cm~ 2 
ijVasisht et al.lll996|) . In comparison, the column den- 
sity fitted to a power-law plus blackbody model for 
XTE J1810-197 was N H = 1.02 x 10 22 cm" 2 . Since 
on physical grounds we strongly prefer the double black- 
body model, which requires Nn — 0.65 x 10 22 cm -2 , 
it seems appropriate to reduce the distance estimate to 
one that is more compatible with the smaller iVu value. 
We now assume that the distance to XTE J1810— 197 
is 2.5 kpc, and explore the consequences of this revi- 
sion. The effect is to reduce the inferred X-ray luminosi- 
ties and blackbody surface areas. For d = 2.5 kpc, the 
surface area emitting the warm blackbody component is 
A\ < 3 x 10 12 cm~ 2 at all times, which is less than 1/6 
the area of a neutron star. This makes it easier to under- 
stand how the pulsed fraction of the softest X-rays can 
be as large as 26%. 

We can also revise the estimate of the total energy 
emitted in the outburst. At a distance of 2.5 kpc, the 
bolometric luminosity quoted in Paper II is £bb2 = 

1.5 X 10 35 d2 5 g-(t-to)/300days ergg g -l w j th respect to 

the initial time to o f the outburst observed by RXTE 
ijlbrahim et alJl2004j) . which corresponds to an extrap- 
olated energy of 3.9 x 10 42 d\ b ergs. Since RXTE is 
only sensitive to the hotter blackbody component, wc 
should add another contribution from the XMM-Newton 
measured warm component, which can be estimated as 

£BB1 = 1-7 X 1034 rf2 5 e-(t-«o)/900 days ergg g -l_ The 

extrapolation required to integrate the contribution of 
the warm component is less certain than for the hot 
component, but it corresponds to an energy of only 
1.3 x 10 42 d\ 5 ergs. The total estimated energy is then 
5.2 x 10 42 d\ 5 ergs, which is comparable to the amount of 
heat assumed t o be deposited in the crust during a dee p 
heating event l)Lvubarskv. Eichler. &: Thompson! 12002). 
or the extra energy stored in an azimuthally twisted ex- 



ternal magnetic field ({Thompson et~alll2002l) . In the fol- 
lowing section, we discuss how the observations might 
probe the actual mechanism of the outburst. 

3.3. Heating or Cooling'? 

Two mechanisms have been invoked to explain the late- 
time "afterglow" of an outburst of a magnetar. The first 
involves surface heating by long-lived currents flowing on 
closed but azimuthally twisted external magnetic field 
lines, and Comptonization of the resulting surface X-ra y 
emission by the same particles ({Thompson et al.l f2002). 
The originating event could be a sudden fracture in the 
crust in response to twisting of the internal magnetic 
field, or a plastic deformation of the crust that gradually 
transfers internal magnetic twist to the external field. It 
is thought that the major outbursts of SGRs are magneti- 
cally trapped fireballs triggered by a sudden fracture, but 
this is not known to apply to XTE J1810— 197 because 
no bursting behavior was seen, and no prior rotational 
ephemeris exists to test for a glitch. In any case, the 
extra energy above that of a pure dipole, A_E, stored in 
a twisted magnetic field external to the star is, 

(4) 

where A0n-s < 1 rad is the azimuthal twist f rom north 
to south hemisphere ( Thompso n et al] 12002'!. This is 
enough to account for the extrapolated 5.2 x 10 42 d\ 5 
ergs radiated if A0n-s is a small fraction of a radian. 

In this model, the X-ray spectrum is a combination 
of surface thermal emission, and Compton scattering 
and cyclotron resonant scattering of the thermal X- 
rays. Compton upscattered flux may result in power- 
law emission that dominates at high energies such as 
that reported above 10 keV for the AXP in Kes 73 
(jKuiper et al. Il2004|) . However, the flat power-law index 
(r ~ 0.94) does not contribute to the softer spectrum, in 
particular below 2 keV. So we consider that the potential 
power-law tail in the fainter XTE J1810— 197 may no be 
detectable yet. The decay timescale is determined by the 
rate at which power is consumed by the magnetosphcric 
currents, which is comparable to that needed to power 
the observed X-ray luminosities of AXPs in general and 
XTE J1810-197in outburst. Furthermore, it is the evo- 
lution of surface heating that basically determines the 
instantaneous X-ray luminosity and its decay, but more 
specific predictions of decay curves have not been made. 

In the second picture, a deep crustal heating event, 
heat deposited suddenly by rearrangement of the mag- 
netic field is gradually transferred by diffusion and ra- 
diation, and the radiated X-rays are purely thermal. 
While the originating event might have been a fracture 
of the crust as in the first scenario, an assumption that 
the heat is evenly distributed throughout the crust re- 
sults in a particular de cay curve from one-dim ensional 
models of heat transfer l)Lvubarskv et al.ll2002f) . In this 
model, the heating is virtually instantaneous, within 
10 4 s, while it is the rate of cond uction cooling that deter- 
mines the observed decay curve. iLvubarskv et al.l l)2002f) 
assumed a deposition of ~ 10 25 ergs cm -3 to a depth of 
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500 m, which is < 1% of the magnetic energy density 
at the surface. If occurring over the entire crust, this is 
~ 1 x 10 43 ergs, comparable to the total X-ray energy 
emitted during t h e dec ay of XTE J1810— 197. However, 
iLvubarskv et alJ l|2002f l showed that the resulting cool- 
ing luminosity would follow a t~ ' 7 power law, which is 
rather slow compared to the observed exponential decay 
of the hot component of XTE J1810-197, while 80% of 
the energy should be conducted into the neutron star 
core rather than radiated from the surface. 

It is not obvious if the behavior of XTE J1810-197 
supports or excludes either of these models. One in- 
teresting possibility is that the slowly decaying warm 
component is the radiation from a deep heating event 
that affected a large fraction of the crust, while the hot 
component of the spectrum is powered by external sur- 
face heating at the foot-points of twisted magnetic field 
lines by magnetospheric currents that are decaying more 
rapidly. It is not yet clear if the flux decay of the warm 
component is primarily an effect of decreasing tempera- 
ture or decreasing area, because the decay has so far been 
modest. It is possible that the decay of the warm com- 
ponent is consistent with a t -0 - 7 power law rather than 
an exponential of n = 900 days. Fits to Lbbi as a func- 
tion of time allows a power-law decay index in the range 
—0.3 to —0.7. In contrast, because of its more rapid de- 
cay, it is evident that the hotter component is declining 
exponentially, in area rather than in temperature. This 
shrinking in area might be easier to understand as a de- 
cay of the currents or a rearrangement of the magnetic 
field lines that are channeling the heating on the surface. 
But there is no good evidence yet of an inverse Compton 
scattered component from the particles responsible for 
that current, at least not at energies < 8 keV. Perhaps 
the solid angle subtended by the twisted field lines in the 
magnetosphere is too small to scatter most of the ther- 
mal photons. Each of the models is challenged in some 
aspect by the observations, but each may still have some 
applicability. Observing the decay through to quiescence 
should help to clarify the situation. 

3.4. Constraints on Emission and Viewing Geometry 

The pulsed light curves offer an additional diagnostic 
of the emission and viewing geometry. If we adopt the 
hypothesis that virtually all of the 0.5 — 8 keV flux is 
surface thermal emission, then the symmetry of the light 
curves and their strict pulse-phase alignment as a func- 
tion of energy argue for a concentric geometry in which 
a small hot spot is surrounded by a larger, warmer re- 
gion. The large observed pulsed fractions are achievable 
in a realistically modeled NS atmosphere that accounts 
for the different opacities of the normal modes of polar- 
ization in a strong magnetic field (jQzelll200ll lOzel et all 
l200l . 

Many geometries were modeled bv lOzel i|200T|) . but 
only total pulsed fractions were reported rather than de- 
tailed light curves. It is important to model the pulse 
profiles completely, since the number of peaks and their 
phase relationship depend in detail on the surface and 



viewing geometry. In addition, the relative effects of fan 
versus pencil beaming, as regulated by the anisotropic 
opacity in a strong magnetic field, allow the number of 
peaks in the light curve to differ from the number of hot 
spots on the surface. Whereas the pulsed fraction of the 
hot component in XTE J1810-197 is 53%, this is not 
high enough to require a single spot. It is possible that 
the flat interpulse region of the triangular, hard X-ray 
pulse is actually emission from an antipodal spot at large 
viewing ang le. lOzel et alJ i|2001l) showed that for large 
magnetic fields, most of the flux would be in a broad 
fan beam, while iThompson et alJ (j2002f) argue that cy- 
clotron resonant scattering in the magnetosphere can sig- 
nificantly reduce this effect. Detailed physical modeling 
of the pulse profiles of XTE J1810— 197 is needed before 
definite conclusions about the emission and viewing ge- 
ometries can be drawn. 

4. CONCLUSIONS AND FUTURE WORK 

Four sets of XMM-Newton observations of 
XTE J1810— 197 spanning a year, clarify several 
new behaviors that might apply to other AXPs and 
AXP-like objects: 

• a two-temperature blackbody model is arguably 
preferred over the standard blackbody+power law 
model, 

• the components of the compound spectrum are de- 
caying at rates that differ significantly, and 

• the component temperatures of the blackbody 
emission are nearly constant, implying that the 
area of emission is steadily decreasing. 

Several outstanding questions remain that can 
be addressed by continuing X-ray observations of 
XTE J1810-197 over the next few years. Will the cur- 
rent spectrum evolve into the prior quiescent one? Is the 
cool quiescent spectrum still present so that it will begin 
to dominate over the warm temperature component? To 
what degree is the quiescent emission pulsed? Why do 
the temperatures hardly change with time, requiring the 
apparent decrease in area? What is the emission geom- 
etry on the NS and what can we learn from it? Future 
observations will test our predictions for the flux decay 
and pulsed fractions and allow a better understanding 
of the relationship between the spectral components and 
the pulse profile. Ultimately, the goal of this study is 
to disentangle the processes of thermal and non-thermal 
emission to which magnetars convert their energy, during 
outburst and quiescence. 

We thank Fred Jansen and Norbert Schartel for pro- 
viding the three XMM-Newton Target of Opportunity 
observations of XTE J1810-197. XMM-Newton is an 
ESA science mission with instruments and contributions 
directly funded by ESA Member States and NASA. This 
research is supported by NASA grants NNG05GC43G. 
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